Quantum spin glass in anisotropic dipolar systems 
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Abstract. The spin-glass phase in the LiHo x Yi_ x F4 compound is considered. At 
zero transverse field this system is well described by the classical Ising model. At finite 
transverse field deviations from the transverse field quantum Ising model are significant, 
and one must take properly into account the hyperfine interactions, the off-diagonal 
terms in the dipolar interactions, and details of the full J = 8 spin Hamiltonian to 
obtain the correct physical picture. In particular, the system is not a spin glass at 
finite transverse fields and does not show quantum criticality. 



1. Introduction 

The study of spin glasses, and in particular the classical transition between the spin- 
glass (SG) to the paramagnetic (PM) phase, have been thoroughly studied since the 
70's p. As interest in quantum phase transitions has grown, the understanding of 
quantum criticality at the SG to PM phase transition has drawn much theoretical 
interest El G] Experimentally, the quantum phase transition is not accessible in 
most spin glasses, where the strong exchange interaction blocks quantum fluctuations 
at available magnetic fields. However, anisotropic magnetic dipolar systems, notably the 
LiHo x Yi_ x F4 compound, have very weak exchange interactions, and seem ideal for the 
observation of such a transition In this system the dipolar interaction and single- 
ion anisotropy terms have magnitude ~ IK and ~ 10K respectively; with spin J = 8, 
appreciable quantum fluctuations are expected already at transverse fields H± ~ IT. 
The strong easy axis anisotropy means that for T <C 10K, the Ho spins truncate to an 
Ising-like doublet; this anisotropy also strongly suppresses all but the longitudinal terms 
of the dipolar interaction. The LiHoF 4 compound can be diluted by exchanging Ho with 
the nonmagnetic Y atom, resulting in the LiHo x Yi_ x F4 compound with any desired x 
j^]. The angular dependence of the dipolar interaction, in the presence of quenched 
randomness, results in frustration. Thus, the system which is ferromagnetic at x= 1 
jH], turns into a SG for x= 0.167, with T c = 0.13K [71 |S]. For the latter dilution, the 
SG to PM transition was studied as a function of T and i/ijTJIEj- Both the linear and 
nonlinear susceptibility were measured. Despite the importance of this experiment, some 
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of its very interesting features long remained without a proper explanation. Questions 
included (i) the reduction of the cusp in the nonlinear susceptibility with decreasing T; 
(ii) the observation that it is much easier to disorder the SG thermally than quantum 
mechanically; (iii) the sharpness of the low T crossover between a PM response and 
slow relaxation; and (iv) the smallness of the low T critical exponent, as observed in 
the nonlinear susceptibility near the transition. 

Theoretically, the LiHo x Y!_ x F4 system was considered to be a good realization of 
the transverse field Ising model (TFIM) in the electronic degrees of freedom. Recently 
EH it was shown that the real system differs from the above model in two 
significant ways, which affect the physics considerably; (a) the hyperfine (hf ) interaction 
between the Ho electronic and nuclear spins is strong, and for x^C 1 in general, and in 
the SG experiments jHj in particular, dominate the physics E] (b) the off-diagonal 
dipolar terms, although effectively reduced, become essential at any finite H± as they 
reduce the symmetry of the model ^Uj • in this paper we show how these features answer 
the first two of the four questions noted above; we also show that a proper understanding 
of the off-diagonal dipolar terms requires going outside the simple TFIM. 

2. Hyperfine interactions 

The Hamiltonian describing the LiHo x Y!_ x F 4 system is given by a sum of crystal field 
Zeeman, hf, and inter- Ho interaction terms: 

H = H C{ + H Z + H M + H iQt . (1) 

The Ho ion has a J = 8 angular momentum. H c f splits the 17-fold degeneracy, leaving 
3 relevant low energy levels - an Ising-like doublet, denoted ||), ||), and a first excited 
state approximately 10K higher in energy, denoted | r z 2 ) [HE]- H z = — J2i 9j^bH-J^ is 
the Zeeman energy, and H int = — J"jf is dominated by the dipolar interaction 

[TB*] . We denote the easy axes by z, and consider -£/j_||x. It is common to neglect all 
but the longitudinal terms in the dipolar interaction, and drop the hf interaction, thus 
obtaining a low energy TFIM effective Hamiltonian: 

H = - E y ii T t T i - A o(^) E *f , (2) 

hi i 

where fj is a Pauli vector describing the two-level effective electronic spin at spatial 
position r = r,, and Ao oc H± 2 /Qq for small H±. However, both the hf and off- 
diagonal dipolar interactions are of crucial importance for the LiHo x Y!_ x F 4 system. 
The Ho atom is a pure isotope I = 7/2 nuclear spin with contact hf interaction 
Hm = AjJ^ili " J%- Due to the strong anisotropy we consider first the longitudinal 
part of the hf interaction = AjI z J z . This term splits each of the states ||), ||) 
into an eightfold multiplet of nearly equidistant levels, with separation ~ 205mK [T2*] 
between adjacent levels (Fig. 1). This splitting, larger than the typical dipolar energy 
and the relevant experimental temperatures jH] , influences significantly the physics of the 
system. The Ising doublet states have now a definite nuclear spin, i.e. I z = —7/2 for the 
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Figure 1. Splitting of the electronic low energy doublet (| and J.) due to the 
longitudinal hf interaction. The GS doublet, a and a have a definite and opposite 
nuclear spin, ±7/2. Transverse magnetic field couples states with the same nuclear 
spin, as is shown by the dashed lines. 



electronic state |f) and I z = 7/2 for the electronic state ||). A transverse field H± then 
couples a = |T, -7/2) with b = ||, -7/2) and a = ||, 7/2) with 6 = |T, 7/2) (see Fig. 1). 
Thus, the transverse field by itself does not induce quantum fluctuations between the 
relevant Ising doublet GSs, but only re-normalizes their effective spin. The transverse 
hf interactions = Aj[I z J z + (J + J~ + I~J + )/2] allow simultaneous changes in I z 
and J z , and thence quantum fluctuations between the Ising doublet states. However, 
as is detailed in Refs. [HUH] this mechanism is negligible for H± <C fl /fiB- Thus, the 
inclusion of the hf interactions results in three energy scales which dictate the behavior 
of the system, and in particular the position of the crossover between the SG and PM 
phases as function of T and H± At zero field, T c (0) is dictated by the dipolar 
interaction Vo, since the nature of the Ising doublet is not important for the classical 
transition. The behaviour of T c (Hj_) with H± is governed by the strength of the hf 
interaction and the crossover at T = is governed by the anisotropy energy Qq, as 
quantum fluctuations between the relevant Ising states become significant only when 
the state | T l 2 ) becomes appreciably hybridized with |T),|I)- All three energy scales 
become apparent in the position of the line separating the SG and PM phases, and in 
particular the relation Vq -C Qq explains the fact that it is much easier to disorder the 
SG thermally rather than quantum mechanically |9| lllj. 

3. Off-diagonal terms of the dipolar interaction 

The longitudinal term of the dipolar interactions, oc JfJ?, has direct matrix elements 
within the low energy Ising states |T, —7/2), ||, 7/2). All other interaction terms involve 
the state | Y l 2 ) in second or higher perturbation expansion. Since flo 3> Vo, one is 
tempted to neglect all but the longitudinal interaction. However, the off-diagonal terms, 
and in particular the terms oc J* JJ become important at H± ^ 0, as they change the 
symmetry of the system. The J z — >• —J z symmetry, while maintained by H± without the 
off-diagonal terms, is destroyed. This reduction of symmetry results in a generation of an 
effective random field at each site [HI]. As a result, within the droplet picture of Fisher 
and Huse and using an Imry-Ma like argument ^3] one obtains jTUj a magnetic 
field dependent finite correlation length. The system with transverse magnetic field and 
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off-diagonal dipolar interactions becomes equivalent to the random field Ising model, i.e. 
a SG in the presence of random longitudinal magnetic field. Thus, the scaling theory of 
Fisher and Huse [Hj predicts an instability of the SG phase to finite H± in our case, 
equivalent to its prediction of the absence of de-Almeida Thouless line in the RFIM. 
Interestingly, to obtain the correct physical picture of the system at finite field one has 
to consider the large Ho spin, going beyond the simplified Ising picture. The reason 
is that it is the fluctuations between each of the single Ho GSs and its first excited 
states, and not the much smaller quantum fluctuations between the two Ising states a, a 
|lUj . that govern the magnitude of the effective random field and the reduction of the 
correlation length. 

Within the scaling picture of Fisher and Huse ^1] the SG at zero transverse field 
has two time-reversed GSs, denoted ip and ip. Each Ho ion in the GS ip is in either states 
a or a, and in the opposite state in the GS ip. Consider a single Ho ion at H± ^ 0. For 
small magnetic fields the fluctuations between the Ising states a and a are negligible 
[Oj. Yet, the energy of each of the Ising states is reduced by an energy proportional 
to H± 2 /Q , due to fluctuations to the relevant excited states (| T^, — 7/2) for a and 
| r z 2 , 7/2) for a). If we choose an arbitrary region in state ip, the energy reduction due 
to the field is just the sum over all spins of the single spin energy gain. With J z — > —J z 
symmetry, the same energy reduction occurs for the state ip. Now consider the effect 
of the off-diagonal dipolar terms (in particular, the term V-^JfJj) on the different 
domains. In second order perturbation, the domain energy shift is given by JU] 

^- — ~n • (3) 

The dipolar terms have randomly the same or the opposite sign to that of the magnetic 
field, and typically, by flipping a domain of N spins one gains an energy of ^U] 



m = c Mv VN 

where j = max J 2 . Comparing this energy gain to the energy cost of flipping the 
domain ~ j 2 VoL 0d where L is the domain linear size, one finds a finite correlation 
length at any H± given bvjlOj 

Importantly, this correlation length depends only on H± and Q Q . In the experiment [S], 
as T is decreased the crossover to the PM phase occurs at higher H±, dictating a smaller 
correlation length and a reduced cusp in the nonlinear susceptibility jTHj. Note that the 
cusp is further reduced due to the renormalization of the effective spin jH] . Interestingly, 
the finite correlation length results in an enhanced transverse field in the x direction 
which was anticipated by the comparison of the experimental and theoretical positions 
of the crossover line between the SG and PM phases jHj. 
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In principle, one could also get a finite correlation length in a spin-half Ising model, 
by introducing a longitudinal interaction oc Vq and a reduced off-diagonal interaction 
oc fttVn|l()j. Similar considerations to the ones above lead to a correlation length 
£ m [V / (a/jBHx)} 1 ^ 3 /^ 6 ^ [TH], coming from quantum fluctuations between the Ising 
doublet states, with a £ depending on Vq. A notable difference between this result and 
the correct one is the size £ at the crossover to the PM phase at T = 0. This crossover 
actually occurs at H± ~ fio for the case of large spin [HIE], leading to £ « 1. However 
in the effective spin-half model the crossover occurs at H± rs Vo and £ ~ 1/q at the 
crossover, emphasizing the inadequacy of the spin-^ model. Importantly, our results 
here are easily generalized to any anisotropic systems, as long as a dipolar interaction 
exists [TUj . 



4. Numerical results 



The scaling relation (J3J) has been checked against Lanczos exact diagonalization (ED) 
computations on finite size clusters ^2]. In order to get closer to the experiment, we 
randomly distribute N moments at the rare earth sites of three dimensional LiHo x Y 1 _ x F4 
diluted lattices, and focus on s = 1 particles with an on-site anisotropy Qq ~ 10K which 
accounts for the crystal field. Therefore, the following transverse field dipolar spin-1 
Hamiltonian 

1 



E 



(6) 



has been diagonalized on LiHo x Y!_ x F4 lattices with x = 18.75% for various sizes, and 
over 10,000 independent random samples for each size. In the perturbative regime, we 
have computed the finite size gap 5E for each sample and the \/N scaling stated in Eq. 
(J1J) is clearly demonstrated, as shown in Fig. [21 (a). 
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Figure 2. Distributions of the finite size gaps 6E rescaled by \f~N and plotted in a 
semi-log scale. Lanczos ED data collected over 10,000 random diluted LiHo x Yi_ x F4 
samples for each size N. Insets: Linear dependence of the disorder average gap (6E) 
vs v^V- (a) Results obtained for the spin-1 Hamiltonian |JBJ with Qo/l^BHt = 100 for 
three different sizes, with x = 18.75%. (b) Results obtained for the spin-i model 
with a = 0.1, jtte-fff = 0.0025 and five different sizes with a dilution x= 1/12. 
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A similar calculation was done for the spin-^ Ising model 



Hi 

2 



E 



\rzz QZ QZ 

2 V v ^i^j 



aV%*SfSj 



(7) 



Numerically speaking, the s = 1/2 problem is easier since it leads to a smaller Hilbert 
space dimension and allows us to check the \/N scaling relation over a broader range of 
sample sizes, as shown in Fig. El(b). However, here too the essential difference between 
the spin-| and spin-1 models is apparent. In the spin-| model, since the significant 
fluctuations are between the Ising doublet states, the strength of the dipolar interaction 
replaces f2 in the denominator of Eq. which results in better convergence of the 
numerics at the smallest sizes. 



5. Conclusion 



Of the four experimental puzzles mentioned in the introduction, the first two are 
explained by our analysis. Since we believe that the Hamiltonian JI} includes all the 
physics relevant to the experiment we expect that the solution to the latter two 
puzzles lies within the same framework of considerations above. 
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